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ABSTRACT 

We use high resolution cosmological hydrodynamical simulations to demonstrate that cold flow gas 
accretion, particularly along filaments, modifies the standard picture of gas accretion and cooling onto 
galaxy disks. In the standard picture, all gas is initially heated to the virial temperature of the galaxy 
as it enters the virial radius. Low mass galaxies are instead dominated by accretion of gas that stays 
well below the virial temperature, and even when a hot halo is able to develop in more massive galaxies 
there exist dense filaments that penetrate inside of the virial radius and deliver cold gas to the central 
galaxy. For galaxies up to ~L*, this cold accretion gas is responsible for the star formation in the 
disk at all times to the present. Even for galaxies at higher masses, cold fiows dominate the growth of 
the disk at early times. Within this modified picture, galaxies are able to accrete a large mass of cold 
gas, with lower initial gas temperatures leading to shorter cooling times to reach the disk. Although 
star formation in the disk is mitigated by supernovae feedback, the short cooling times allow for the 
growth of stellar disks at higher rcdshifts than predicted by the standard model. 
Subject headings: galaxies: evolution — galaxies: formation — methods: N-Body simulations 



1. INTRODUCTION 

The classic model for galaxy formation within the cold 
dark matter (CDM) paradigm assumes that gas within 
galaxy halos initially reaches the temperature of the viri- 
alized halo. During the collapse of the halo, the pro- 
cess of violent relaxation of the dark matter shocks the 
gas component to the virial temperature. Any subse- 
quently accreted gas is also shock ed to the virial temper- 
ature as it enters the viria l radius (Rees fc Ostrikei]fl977l : 
ISiig 119771 : IWhite fc Reei lT9781. This assumption has 
been routinely adopted by many analytical works (e.g., 



Whit.e fc Frenklll991t iKauSmann ei^LlfTool [CoIeeTaLl 

fl998l: 



19941: lAvila-Reese et al.l Il998t ISomerville fc Primac 
1999"; 'Stringe r fc BensonI l200l iKampakoglou fc Sil 
2007; Silk 200f~ Given a density profile, a cooling ra- 
dius can be calculated, inside of which gas can radi- 
ate away its energy while conserving its specific angular 
momentum to form a centrifug ally supported disk that 
grows from the inside out (e.g., [ Fall fc Efstathiou" 19801; 
Dalcanton ct a l1ll997t iMo et al.i 1998: Mao ct al. 1998; 
van don Bosch 2GG0l: lHatton et al.ll2003t lSomcrvill e et all 
2008). In general, these semi-analytic models (SAMs) 
have been successful in reproducing many of the present 
day properties of galaxies. 

Some studies have called this virial temperature as- 
sumption into question, however. Recent theoreti- 
cal work has found that, in fact, the majority of 
baryons may never reach t he virial temperat u re as 
they accrete onto galaxie s (iBirnboim fc Dekell l2003t 
Dgkcl fc Birnboim 2006; K eres et al.ll2005l : lOcvirk et al I 
2008t IKeres et al.ll2008t iDekel et al.ll2008f ). iKeres et al.l 
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(|2005| ): IKeres et al.l (|2008l hereafter K05 and K08, re- 
spectively) examined how gas gets into galaxies using 
a N-Body -I- Smoothed Particle Hydrodynamics (SPH) 
simulation of a cosmological volume. They found that 
galaxies with halo masses below a few 10^ M© are dom- 
inated by "cold" gas accretion. This cold gas is never 
shock heated to the virial temperature of the galaxy, and 
is thought to be accret ed predominantly a long filaments 
in the cosmic web (e.g.. lMassev et al.ll2007D . Further sup- 
port f or this new paradigm has come from lOcvirk et al.l 
()2008l hereafter OPT08) using a similarly large cosmo- 
logical volume simulated with an adaptive mesh refine- 
ment (AMR) code. 

For galaxy halos to shock heat infalling gas, the cool- 
ing rate of gas behind the shock must be slower than the 
compression rate of the infallin g gas (jBirnboim fc Dekell 
120031 : IDekel fc Birnboiml[200l l. This requirement leads 
to two regimes in which gas is able to avoid shock heat- 
ing, and instead cools onto the galaxy on a free fall 
timescale instead. First, galaxies below a certain crit- 
ical mass are not massive enough to support a stable 
shock because their dynamical time is shorter than the 
cooling time of the halo gas. OPT08 find that the 
halo mass where cold and shocked accretion contribute 
equal ly is in good agreement \y ith the theoretical analy- 
sis of IDekel fc BirnboimI ()2006l hereafter DB06), with a 
critical mass of IQ-'^^-^ M0. This critical mass is roughly 
independent of redshift. 

Second, even when a shock is present, cold gas ac- 
cretion can occur along filaments that penetrate deep 
inside the hot halo. The ability of filaments to 
feed galaxy growth has been noted for some time 
(!Katz et al.''l993'; 'Katz & White' ■1993"; 'Katz ct al. 199 
Bertschinger fc Jain 1994; Bond ct al. 1996; Shen ct al^ 
boOetlHarford et al.ll2008( ). K05 find that their cold flow 
gas accretion is largely along filaments. These filaments 
are efficient at funneling gas to the galaxy, and can bring 
material to the galaxy from a much larger radius than 
predicted in the standard spherical accretion model. The 
mass contributed by filaments is a strong function of red- 
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shift (DB06, OPT08). 

Note that in the first case above, when a galaxy is 
in the low mass regime incapable of supporting a sta- 
ble shock, cold gas accretion can occur both in a quasi- 
spherical fashion and along filaments, while in the sec- 
ond, high mass regime in which a hot halo exists, cold 
gas accretion occurs primarily along preferred directions 
in filaments. 

Galaxies above the critical mass for shocking appear 
to be dominated by hot gas accretion as in the standard 
model, but they may have been fed by cold flow filaments 
even after a shock develops, particularly at high z, and 
they are also built in part by the mergers of smaller galax- 
ies that were dominated by cold flows. Thus, cold flow 
gas accretion is an important component of galaxy for- 
mation that has remained largely ignored until recently. 

Previous work on cold flow gas accretion in simula- 
tions have addressed the issue using large cosmological 
volumes, furnishing a statistical description of how galax- 
ies get their gas. This provides a useful starting point 
for modifying semi-analytic models of galaxy formation. 
However, with large volumes and many simulated galax- 
ies, numerical resolution for any individual galaxy must 
be sacrificed. Hence, these simulations were unable to re- 
solve the internal structure of galaxies and investigate the 
impact of cold flows on the formation of galactic disks. 
OPT08 follow gas as it enters the virial radius of galax- 
ies, and examined the temperature history of that gas 
as it flows toward the central galaxy, down to 0.1 Kvir- 
However, limited resolution prevented them, or K05 or 
K08, from being able to resolve disks in their galaxies. 
Thus, the impact of cold flows to disk assembly, and in 
particular the growth of the stellar disk component, has 
yet to be addressed. 

Rather than sample a cosmological volume at low res- 
olution, in this paper we investigate gas accretion in sev- 
eral very high resolution simulations, in order to study 
the impact of cold gas accretion on the assembly of stellar 
disks. Each of these galaxies has at least a million par- 
ticles within the virial radius by redshift zero. This high 
resolution provides two advantages over previous simula- 
tion work on this subject. First, it allows us to resolve the 
internal structure of our galaxies and their disks. Second, 
we overcome artificial broadening of shocks in SPH so 
that we are able to identify entropy increases associated 
with accretion shocks, allowing us to explicitly search for 
shocks rather than simply using a temperature threshold 
as in previous works (K05, K08, OPT08). 

We are thus capable of examining the temperature and 
entropy history of accreted gas that eventually settles 
to the disk and forms stars. These simulated galaxies 
are run within a ACDM context and span two orders of 
magnitude in galaxy halo mass, allowing us to study the 
role of cold fiows in galaxies from 0.01 L* to L*. 

This paper is outlined as follows: In Section [2] we 
briefly discuss the simulations used in this study. In Sec- 
tion [3] we describe our method to identify shocked and 
unshocked gas accretion to galaxies, and flnd that the 
fraction of accreted gas that is shocked is a strong func- 
tion of galaxy mass, in agreement with previous works. 
In Section [4] we extend the results of Section [3] to in- 
vestigate the role that unshocked, cold gas accretion has 
in building the stellar disk of galaxies, as a function of 
galaxy mass. We briefly compare our results on disk 



growth to analytic models of disk galaxy formation. In 
Section [5] we discuss the various effects of our numerical 
scheme on the results, and conclude that our results are 
robust (those readers interested strictly in the science re- 
sults may wish to skip this section). Finally, in Section 
[6] we discuss the implications of our results for existing 
analytic models of disk galaxy formation, and conclude. 



2. THE SIMULATIONS 

The simulations used in this study are the culmina- 
tion of an effort to create realisti c disk galaxies. As dis- 
cussed in iGovernato et al.l ()2007l ) , at z = disk galax- 
ies in our simulations overcome the angular momentum 
catastro phe and resolve the in ternal structure of disks 
(see also IGovernato et al.ll2008L which shows that galax- 
ies in this paper lie on the TuUy-Fisher relation). This 
crucial improvement allows us to undertake the current 
investigation, and is made possible both by dramati- 
cally increased numerical resolution and by the inclu- 
sion of a simple but physically motivated recipe to de- 
scribe star formation (SF) and the effects of subsequent 
super nova (SN) feedback (full details in IStins on et al.l 
I2006f) . Our adopted SN feedback and cosmic UV back- 
ground (Haardt & Madau 1996) drastically reduce the 
number of galaxy satellites containing a significant stel- 
lar population, av oiding the well known "substructure 
problem" (W hite fc Frenklll99l|: iKauffmann et al.|[l99l 
iQuinn et al.lll996t iMoore et al.N l999'). Our SN feedback 
is also effective at regulating SF, preventing gas from ar- 
tificially cooling too early and forming stars, and hence 
we produce smaller bulg es and halos in subsequent merg- 
ers (jBrook et al.ll2004aJ ). This regulation by SN also re- 
produces the observed stellar mass-metallicity relation- 
ship for galaxies, both at the present time and a t high 
redshift (|Brooks et al.ll2007l : iMaiolino et al.ll2008D . At z 
= 3, these simulations reproduce the metallicities, inci- 
dence rate, and column density distributi ons of observed 
damped Lyman alpha absorbers (DLAs iPontzen et al.l 
2008). 

Four simulated disk "field" galaxies were chosen to 
cover two orders of magnitude in halo virial mass, from 
3.4x10^° Mq to 3.3x10^2 j^^^ ^he th ree most massive 
galaxies have been discussed in detail in IGovernato et aTl 
(2007, 200^, i.e., DWFl, MWl, GALl), but they have 
been rerun at eight times the mass resolution and twice 



the spatial resolution for this study. All galaxi es are sim- 
ulated using Gasoline (jWadslev et al.ll2004l ) in a fiat, 



A-dominated cosmology. The three most massive galax- 
ies were run using a WMAP year 1 cosmology, while a 
WMAP year 3 cosmology was adopted for the newest, 
lowest mass galaxy (hereafter H579). The three least 
massive galaxies were originally drawn from uniform vol- 
umes ranging from 25 Mpc to 50 Mpc, while GALl was 
selected from a 100 Mpc box. These galaxies were then 
resimulated at higher resolution using the volume renor- 
malization technique (jKatz &: Whitd[l993f l . which allows 
us to resolve fine structure while capturing the effect of 
large scale torques. We adopt th e same SN and SF ef- 
ficiency parameters as in lGover"riato et al. ( 200^ and a 
Kroupa initial mass function (jKroupa et al.ill993l ). 

Characteristics of each galaxy at z = are given in 
Table 1. Galaxies and their parent halos were identified 
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TABLE 1 

Simulated Galaxy Properties 



simulation 




R-Tjir 


^ circ 


T 


A 




e 


N within R„ir 




M0 


kpc 


km/s 


K 






kpc 


dm+star+gas 


H579 


3.4x10^° 


85 


40 


5.1x10"' 


0.03 


>5 


0.17 


~9.5xl05 


DWFl 


1.4x10" 


135 


70 


1.3x10^ 


0.01 


1.5 


0.15 


-5.3x10*^ 


MWl 


1.1x10^2 


270 


130 


5.1x10^ 


0.05 


2.5 


0.3 


~4.8xl0'^ 


GALl 


3.3x10^2 


385 


190 


1.1x10^ 


0.02 


2.8 


0.3 


-3.7x10*^ 


H277 


7.1x10" 


230 


115 


3.8x10^ 


0.03 


2.5 


0.35 


~2.3xl0'^ 


MWl.lo 


1.1x10^2 


265 


130 


5.0x10^ 


0.07 


2.5 


0.6 


~5.8xl05 


MWl. ad 


1.1x10^2 


270 


130 


5.1x10^ 


0.07 


2.5 


0.6 


~3.0xl0'^ 



using AHF^ (iKnebe et al.l[200lt iGill et al.l[200l . AHF 

determines the virial radius, Hvir, for each halo at each 
output time step using the over density criterion for a flat 
universe following [Grosd (|1997D . AHF identifies all halos 
that exist in the high resolution region of our simulations 
at each output time step, and the individual particles 
that belong to them. We can thus track a single gas 
particle through the entire simulation, as it is accreted 
to a galaxy, as that galaxy later merges with another, 
and as the gas particle forms stars. 

AHF allows us to identify the main galaxy at each out- 
put time step from z = to higher z, and the particles 
that belong to it. In this paper, we identify the progen- 
itor halos of the main galaxy all the way back to z = 
6. Thus, as can be seen in Table 1, particles are traced 
back to a higher z than the last major merger for these 
galaxies. For example, our MWl galaxy undergoes a last 
major merger at z ~2.5, and is easily identified at lower 
redshifts. At higher z, we follow the most massive of the 
t wo progenitors to study the galaxy's history. 

iGovernato et all (|2007[ ) and lReed et"al] (|2003D showed 
that the mass function in their simulations converged at 
a minimum of 64 dark matter (DM) particles per halo. 
Thus, we require a minimum of 64 DM particles for halo 
identification in AHF. Halos with fewer particles are con- 
sidered as "smooth" gas accretion instead, but in practice 
these halos are "dark" and have no bound baryons due to 
heating from our UV background. The UV background 
turns on at z = 6 and by z = 4 has a significant effect 
on the lowest mass halos, heating their gas so that only 
halos with ^100 or more particles are of sufficient mass 
to retain bound gas particles (with velocities below the 
escape velocity of the halo). Therefore, unresolved halos 
with fewer than 64 DM particles should contribute little 
or no bound gas, and have no significant impact on the 
gas accretion study in this paper. Resolution effects are 
explored further in Section [5.41 

Our highest (lowest) resolution runs have particle 
masses of 4.5x10" Mg (2.9x10^ M©), 1.6x10'' Mq 
(4.3x10^ Mq), and 4.7xl03 Mq (l.SxlO'^ Mg) for dark 
matter, gas, and stars, respectively, and a force resolution 
of 0.15-0.3 kpc. These galaxies were simulated to have a 
similar dynamical range and similar number of particles 
within the virial radiu s. For a more detailed discussion 
of the simulations, see IGovernato et al] (|2007f ). 

We examine the temperature, entropy, and accretion 
histories for gas particles accreted to the main galaxy in 
these four simulations. For all cases, each gas particle's 

^ AMIGA'S Halo Finder, available for download at 
[http://www.aip.de/People/aknebe/AMIGAl 



properties are determined at z = 6, z = 5, and z = 4, 
after which the values are sampled approximately every 
320 Myr until z = 0. The gas in these simulations is 
allowed to cool via atomic cooling, adopting the cooling 
function for a primordial (metal free) gas. We show in 
Section [5.5.21 that adding metal line cooling to the sim- 
ulations has no effect on the results presented in this 
paper. 

Major mergers have been defined as having a ^3:1 
mass ratio, or lower. The galaxies in this sample tend 
to have relatively quiescent merging histories, with the 
last major merger redshifts being greater than 1.5. This 
similarity is by design, in order to compare gas accretion 
histories of galaxies with similar pasts. In a future paper, 
we will explore gas accretion and disk assembly in a set 
of L* galaxies (^10^^ at z =0) that span a broader range 
in space in spin parameters and in last major merger 
redshifts. 

In addition to the four main galaxies studied here, a 
few other simulations have been run to investigate nu- 
merical convergence. MWl is explored at half of the 
dynamical resolution (MWl.lo, the same resolutio n sim- 
ulation discussed at length in IGovernato et all [2007). 
This lower resolution run is also studied with cooling 
turned off (the adiabatic case, MWl. ad) and with metal 
line cooling added (MWl.mc). An additional L* galaxy, 
H277, with M„r = 7.1x10" Mq, is used to investigate 
the effects of output time step resolution on our temper- 
ature and entropy criteria. 

3. GAS ACCRETION TO THE GALAXY HALO 

Before investigating the role of shocked versus un- 
shocked gas in our simulations, we first want to identify 
the gas that has been accreted smoothly to our galaxy 
halos, as opposed to gas that has been accreted from 
other galaxy halos. It is this smoothly accreted gas that 
will be available to potentially shock as in the classic 
model. 

3.1. Smooth vs Clumpy Accretion 

We first find all gas particles that have ever been identi- 
fied as part of the main galaxy, back to z = 6. AHF allows 
us to identify those gas particles which have been part 
of other galaxies before merging with the main galaxy. 
Those gas particles that are already in the main halo 
progenitor at z = 6 (i.e., the earliest step at which gas 
particles are traced) are labeled as "early" accretion. We 
adopt an extreme definition and label as "clumpy" accre- 
tion any particles that have at any output step belonged 
to a galaxy halo other than the main galaxy we're con- 
sidering. The majority of this gas will be accreted in 
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merger events (note that in this definition, even those 
particles that come from the smaller progenitor in a last 
major merger get labeled as clumpy; it includes both the 
hot and cold gas accreted in this merger) , though a small 
fraction of this clumpy gas may also include gas that was 
tidally stripped or blown out of a galaxy halo and later 
accreted to the main galaxy we're considering. The re- 
maining gas that does not meet our "clumpy" criterion 
is labeled as "smoothly" accreted gas. 

Fig. [T] shows the fraction of gas accreted since z = 
6 that has been identified as "clumpy" (green, at top 
of bars) or "smooth" (the blue and red bars, combined; 
this separation of smooth is discussed in detail in sec- 
tion [321) 7 i'^ order of increasing galaxy mass for our four 
main galaxies. The total halo mass (DM-|-stars-|-gas) of 
each of the four galaxies is listed below their respective 
bar. 

The key feature of Fig. [T] is the dominance of smoothly 
accreted gas in the formation of these galaxies, despite 
our extreme definition for "clumpy" accretion. The 
first two columns of Table [2] list the ratios of "smooth" 
and "clumpy" gas accretion to total gas accreted to the 
galaxy since z = 6. In all cases, ~75% or more of the gas 
accreted to the galaxy since z = 6 has never belonged to 
another galaxy halo (the adiabatic case is discussed in 
section 15.5.11) . 

Because we resolve halos down to the mass limit where 
the UV background suppresses galaxy formation, the ma- 
jor numerical uncertainty in our smooth gas accretion 
fractions is due to our finite output time interval (320 
Myr). Due to output time step resolution, our fractions 
in Table [2] for smoothly accreted gas will be an upper 
limit. The issue of output step resolution is explored 
below (section 15. ip and it is shown that, for a Milky 
Way type galaxy, increasing time resolution will identify 
a higher fraction of "clumpy" material, but not by more 
than 10%, so smoothly accreted gas still dominates the 
formation of this galaxy. 

The dominance of smooth gas accretion in the baryonic 
growth of g alaxies agree s with previous s im ulation results 
Keres et al. 2005; Mu rah et all I2002D . iMurali et all 



TABLE 2 

Gas Accretion Properties 



(2002) looked at this question, in particular. Their mass 
resolution allowed them to fully examine gas accretion in 
galaxies with masses greater than 5.4x10^° in halo 
mass, and they found that galaxies above this mass pre- 
dominantly gain mass via smooth accretion rather than 
in mergers. In fact, they find that mergers account for 
no more than 25% of mass accretion at redshifts greater 
than 1, and no more than 35% of growth at 2 < 0.5. 
K05 and K08 also concluded that, in a global sense, the 
growth of galaxies is dominated by smooth gas accretion 
rather than mergers. 

Fig. [1] shows that, in general, there is a trend in the 
amount of "clumpy" accretion as a function of galaxy 
mass. This figure suggests that as galaxy masses get 
smaller, the amount of material accreted from other 
galaxy halos decreases. This is at odds with the merger 
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iFakhouri & Ma 


I2008D. How- 



ever. Fig. [T] represents only accreted gas, and in practice 
there is a decrease in the baryon fraction of halos with 
decreasing halo mass. The reionization of the universe 
acts to unbind baryons from halos with Vcirc below '^30 



galaxy 


^smooth 


fclumpy 


fshock 


funshock 


H579 


0.96 


0.04 


0.00 


1.00 


DWFl 


0.87 


0.13 


0.10 


0.90 


MWl 


0.73 


0.27 


0.65 


0.35 


GALl 


0.75 


0.25 


0.79 


0.21 


H277 


0.73 


0.27 


0.44 


0.56 


MWl.lo 


0.78 


0.22 


0.62 


0.38 


MWl.ad 


0.49 


0.50 


0.84 


0.16 



Note. — The first two columns, famooth and 
fclumpy I a^re the fractional quantities with respect 
to all of the gas accreted to the galaxy since z = 
6. By definition, fsmooth and fciumpy sum to 1.0. 
On the other hand, the last two columns, fghock 
and funshocki are the fractional quantities with 
respect to all of the gas ever smoothly accreted 
to the galaxy, and sum to 1.0. 



km/s, and the fraction of unbound baryons increases as 
halo mass decreases further. Our UV background mimics 
this effect of reionization. Gas loss due to SNe feedback 
also becomes increasingly effective at these lower masses. 
The combined effect is for halos to become increasingly 
dark matter dominated toward lower halo mass (e.g., 
Bullock et a l. 2000; Hocft ct al. l2006l: iGovernato et all 
2007t iBrooks et al. ,2007 : Okamoto et al.l I2008D . As 
smaller halos merge with even smaller, more dark mat- 
ter dominated halos, the amount of baryonic material 
accreted from other halos decreases with smaller halo 
masses, as seen in Fig. [T] 

Fig. [2] shows the accretion rate of gas to the virial 
radius as a function of time for each of the four main 
galaxies. Lines denoting smoothly accreted gas and 
"clumpy" accretion gas are shown in purple (solid) and 
red (dashed), respectively. Major mergers are conspicu- 
ous for all but the lowest mass galaxy, II579, for which 
the last major merger occurred at a redshift greater than 
5. The obvious result here, again, is that smooth accre- 
tion of gas that has never belonged to another galaxy 
halo dominates the gas accretion history of all of these 
disk galaxies. Fig. [2] demonstrates that this is true at 
all redshifts, and not just for the cumulative case that is 
shown in Fig. [TJ 



3.2. Shocked vs Unshocked Gas 

Having now identified two subsets of gas accretion to 
the virial radius, "smooth" or "clumpy," we now focus on 
the gas that has been smoothly accreted. In the classic 
analytic picture of disk galaxy formation, this smoothly 
accreted gas will be shocked as it enters the virial ra- 
dius of the galaxy. DB06 showed, instead, that there is 
a critical mass below which galaxies are unable to sup- 
port a stable shock at the virial radius. If a galaxy is 
below this mass, then gas may not shock until it reaches 
the disk, where the density contrast is high. In massive 
galaxies, even after a stable shock develops and propa- 
gates out through the galaxy's halo, filaments of cold gas 
may penetrate inside the shock radius. These filaments 
may shock when they reach a radius closer to the disk. 

3.2.1. Shock Definitions 
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Fig. 1. — This plot shows the fraction of gas that has been accreted at the virial radius since 2 = 6 that has been accreted as "clumpy" 
(green, top), "unshocked" (blue), and "shocked" (red) gas. Togeth er, t he unshocked and shocked gas make up the total of smoothly 
accreted gas that never belonged to another galaxy halo (see Section l3.2l l. The total halo masses (in M©) of each of the four galaxies are 
listed below their respective bar. 



Smoothly accreted gas is followed until it reaches a ra- 
dius of 30 kpc^ from the galaxy center (or Kmr at high 
z, whichever is smaller). We stop following gas particles 
at 30 kpc to avoid contamination from feedback effects 
that occur within our galaxy disks that also lead to a 
strong increase in entropy and temperature, mimicking 
the shocks we wish to identify. Similarly, we exclude 
"clumpy" gas from consideration in the remaining anal- 
ysis, as feedback from merger induced star formation may 
mimic shocks. 

By definition, a particle undergoes an entropy increase 
when it encounters a shock. Previous work that searched 
for shocked gas within simulations used a temperature 
criterion rather than an entropy criterion. Our high res- 
olution simulations allow us instead to search for the en- 
tropy increase associated with shocks. We searched all 
smoothly accreted gas for an entropy jump that corre- 
sponded to a Mach number of 3 or more. However, we 
found that not all of this shocked gas was associated 
with an increase in temperature to near the virial tem- 
perature of the halo (we discuss in Section 15.21 that this 
is unlikely to be a resolution effect). For the purposes of 
this work, we are most interested in comparing to semi- 
analytic models of galaxy formation, which adopt the 
assumption that gas starts at the virial temperature of 
the halo before cooling to the disk. Thus, our adopted 
definitions below include a temperature criterion. 

At each time step, we determine the densities and 

All distances given in this work are physical, unless specified 
as comoving. 



temperatures for each gas particle and compare to the 
properties of the main halo at that step. In the limit 
of a strong shock, it follows (see also DB06) from the 
Rankinc-Hugoniot shock jump conditions for a singular 
isothermal sphere: 

P shock = 4po (1) 

and 

Tshock ^ ^/^^vir: (2) 

where po is the gas density prior to encountering the 
shock, Pshock is the post-shock density, and Tyir is the 
virial temperature of the galaxy halo that the gas en- 
counters. 

We require that a gas particle must reach 3/8 Tyir 
to be labeled as "shocked." In practice, however, due to 
the heating of the intergalactic medium (IGM) by our UV 
background, 3/8 Tyir of low mass halos is below the tem- 
perature floor of the IGM. Alternatively, this means that 
the entropy of these halos is below the entropy floor of 
the IGM (discussed further in the next section) , where we 
adopt the entropy definition of S = logio(T^-^/p). Thus, 
for a particle to be labeled as "shocked," it must reach a 
minimum entropy threshold in addition to a temperature 
(3/8 Tyir) threshold. This entropy criterion requires 



Sshock > logio[{3/8Tyrr) ■ /4po] 



(3) 



For each smoothly accreted gas particle at each time 
step we ask what minimum change in entropy, AS, is re- 
quired for the particle to go from its initial entropy, Sq 
at step t, to Sshock by the next time step, t+St (where. 
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Fig. 2. — This plot shown the gas accretion rates at the virial radius for "smooth" and "clumpy" accretion as a function of time for each 
of the four main galaxies. The solid, purple line shows the mass in gas that has been accreted smoothly, while the dashed, red line shows 
the mass in gas that has been accreted from other galaxy halos. The present epoch, 2 = 0, occurs at 13.7 Gyr. Narrow, rapid increases in 
the accretion rate are characteristic of mergers. 



again, St is ~320Myr at most redshifts). If the entropy 
at the next time step is in fact > Sshock, and the temper- 
ature of the particle at that next time step is > 3/8 Tyir, 
we count that particle as having undergone a shock. 

In summary, to identify shocked gas we require that a 
gas particle undergoes a strong increase in both entropy 
and temperature, using two simple rules: AS > Sshock- 
So, and T shock > 3/8 Tyir- Gas particles which meet 
these criteria are labeled as "shocked," and the remain- 
ing smoothly accreted gas particles are then identified as 
"unshocked." If filaments undergo either strong shocks 
or a series of weak shocks within the halo, our definitions 
will identify both of these processes. 

Finally, we restrict our label of "shock" to those parti- 
cles that shock at the time they enter the main galaxy's 
virial radius, or afterward. We note that many of our 
particles do in fact undergo entropy jumps which meet 
our definition of shocks at earlier times. It is common 
for these earlier shocks to occur while the gas particle is 
still at distances of several Kmr from the main galaxy, 
and at times many Gyr before they are ultimately ac- 
creted. As we presume these shocks have nothing to do 
with galaxy accretion (our focus here), we ignore them 
(and note that these shocks are likely to occur when the 
UV background turns on, or in collapse in filaments). 

Fig. [1] shows the results of our search. The smoothly 
accreted gas has been divided into "shocked" (red) and 
"unshocked" (blue), based on the temperature and en- 
tropy criteria described above. Immediately obvious is 
that there is a strong trend for the amount of shocked 
gas within a galaxy to increase with virial mass. The last 



two columns of Table [5] quantify the fraction of shocked 
and unshocked gas mass to the total smoothly accreted 
gas mass. When considering only the fraction relative 
to smoothly accreted gas (rather than all gas), the trend 
with mass in the table is even stronger than originally 
appears in Fig. [TJ 

3.2.2. A Mimmum Halo Mass for Shocking 

iBirnboim fc Dekell (|200l . K05, BD06, K08, and 
OPT08 found that there is a minimum halo mass re- 
quired to support a stable shock. Similarly, a minimum 
halo mass to shock is built into our definition, which can 
be deduced from Eqn. [31 as follows. 

Gas particles in the intergalactic medium will be at 
the mean baryonic density for that redshift, ph. How- 
ever, their temperature will be set by the ionizing UV 
background radiation.^ Thus, these particles have some 
entropy, Smean ■ If particles with Smean are to shock and 
reach Sshock as they enter the galaxy's virial radius, their 
densities must jump to 4ph, and their temperature must 
be 3/8 Tvir- For low halo masses, though, with low Tvin 
Sshock is less than Smean- If this is the case, our shock 
definitions would lead to the false conclusion in which 
every accreted gas particle appears to be shocked. Effec- 
tively, the sound speed in the infalling gas is comparable 
to the infall velocity so there is no opportunity for a 
supersonic-subsonic transition. Clearly, we must require 

* The temperature of our UV background was quantified by 
running a simple gas and dark matter run with no star formation, 
no gravity, and cooling on. 
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Sshock > Smean for gas to be physically shocked. 

This requirement then leads to a minimum halo mass 
that is capable of shocking, as a function of redshift. 
For our lowest mass galaxy, H579, Sshock never becomes 
greater than Smean- Hence, all of the smoothly accreted 
gas in H579 is identified as unshocked gas. DWFl does 
not shock gas until it reaches a halo mass of 1.1x10^^ 
Mq at z — 0.6. This late shock development is seen in 
detail in Fig. [Sj which is the same as Fig. [21 but now 
with smoothly accreted gas divided into shocked and un- 
shocked gas. Meanwhile, our two most massive galaxies 
(MWl and GALl) reach large enough masses to shock 
very early on. For MWl, this happens at z = 4 when the 
galaxy is 8.9x10^° Mq. GALl is 1.6x10" M© at z = 6 
(as far back as we trace gas), at which time it is already 
capable of supporting shocks. 

Although these four galaxies do not provide a statisti- 
cal sample, the mass at which their shocks exist are in 
good agreement with previous results. For DWFl and 
MWl, the masses listed in the previous paragraph are 
the masses at which a shock is capable of developing, and 
thus when a hot halo might develop. GALl is already 
capable of shocking at z = 6, so this mass is an up- 
per limit to when a stable shock develops. OPT08 show 
that shocking near the disk radius begins to occur when 
dark matter halos reach masses of a fewxlO^*^ Mq to 
10^^ MQ(see their figure 4), in good agreement with the 
masses found for our galaxies. We note that the shock 
definition of these authors differs from ours, but find be- 
low (Section 15. 3p that our two definitions yield similar 
results. 

These four galaxies were chosen to have their last ma- 
jor mergers at z > 1.5. If we had instead examined 
galaxies with late major mergers, they may not cross the 
shock threshold until low redshifts, and have even higher 
fractions of unshocked gas accretion. Our galaxies here 
assemble the majority of their mass early, and then qui- 
etly accrete gas afterward, meaning that they effectively 
cross the mass threshold for stable shocks early. Thus, 
by adopting galaxies with similar merger histories, we 
guarantee that we are comparing galaxies across a range 
in mass, but with other properties being similar. 

Because of their range in final mass, each galaxy crosses 
the shock threshold at differing times. Galaxies that 
cross this threshold sooner, i.e., more massive galaxies, 
therefore have a longer period of time to accrete shocked 
gas. This extended time capable of shocking is partly 
responsible for the well defined trend of accreted shock 
fraction with increasing mass, in that these galaxies have 
similar histories and thus display this trend sharply. 
However, larger volume cosmological simulations that are 
capable of statistically sampling a large number of galax- 
ies also show this trend (K05, K08, OPT08). 

3.2.3. The Temperature History of Shocked and Unshocked 

Gas 

Fig [4] shows histograms of the maximum temperature 
attained by all the gas that is smoothly accreted by our 
four galaxies. The temperature shown here is the max- 
imum temperature the particle reaches before we stop 
tracing it (i.e., until the particle reaches 30 kpc of the 
galaxy's center). The solid line shows all of the smoothly 
accreted gas, while the dotted lines show the component 
of that smooth accretion that has remain unshocked, and 



the dashed lines show the shocked component. 

In general, the shocked/unshocked gas is well divided 
into hot/cold as well. H579 has no shocked gas, but 
for our two most massive galaxies, MWl and GALl, 
the two components are well separated in temperature. 
This separation is not altogether surprising, as part of 
our requirement for being identified as shocked is that 
the particle reach a minimum temperature (3/8 T„ir) as 
well as have a sudden entropy increase. The exception is 
DWFl. While the unshocked gas extends to much lower 
temperatures than the hot gas, it overlaps the shocked 
gas. This result is due to the requirement mentioned 
above that the shocked entropy of the galaxy must be 
above the ambient, intergalactic entropy (Smean, set by 
the mean baryon density and the UV temperature floor). 
DWFl does not cross this threshold until z ^0.6. Conse- 
quently, there exists gas that reaches temperatures above 
3/8 Tyir of DWFl at higher redshifts, and is not included 
as shocked. 

It is difficult to say immediately whether this gas 
should be counted as shocked gas. Note that even the 
two most massive galaxies, which are capable of sup- 
porting a shock over basically their entire traced histo- 
ries, include a portion of gas that reached high temper- 
atures (above a few 10^ K) yet is labeled as unshocked 
because it does not undergo a strong entropy increase. 
Previous works on this subject have attempted to iden- 
tify shocked and unshocked gas by using only a temper- 
ature cut (e.g., K05 and OPT08 used a cut of 2.5 xlO^ 
K). In that case, our hot, unshocked gas would have been 
included as shocked gas. We will investigate this alter- 
native definition to search for hot and cold flows later in 
the discussion (Section [5?3|) . 

Despite the small inclusion of hot gas in the unshocked 
definition, we find that our definitions do an excellent 
job of separating gas spatially, into gas accreted along 
preferred directions in filaments (unshocked) and gas ac- 
creted from a more spherical distribution (shocked) . This 
distinction was also found by K05 and OPT08, using 
their temperature definition. Fig. [5] demonstrates this 
for our MWl galaxy. At z = 4 and 3, the left hand col- 
umn shows in green the gas particles that will eventually 
be (or have already been) accreted to the galaxy but re- 
main unshocked, while the right hand column shows in 
green those gas particles that will be accreted (or already 
have been) to the galaxy and shocked. The underlying 
colors are simply a gas density map, for reference (black 
is least dense, white is densest structure). Each frame is 
~1100 comoving kpc on a side. Shown as a faint green 
box is the virial radius of our MWl galaxy at each time. 
At z = 4, the galaxy has only just crossed the entropy 
threshold to allow shocked gas, so that it is only accret- 
ing unshocked gas at this step. However, by z = 3 the 
galaxy is capable of accreting shocked gas, but is still 
accreting a substantial amount of unshocked gas in fila- 
ments penetrating deep inside of the virial radius. This 
demonstrates that the two modes of accretion can occur 
simultaneously at high z. 

The immediately obvious feature in Fig. [S] is the clear 
distinction in spatial distribution at these high redshifts, 
before most of the gas has been accreted. The un- 
shocked gas is collapsing into filaments that are feeding 
the galaxy's growth. Because these filaments are capable 
of penetrating deep inside the virial radius, this cold gas 
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Fig. 3. — This plot is the same as Fig. |2]except that the "smooth" gas accretion rate has been divided further into "shocked" (red, solid 
line) and "unshocked" (blue, dot-dashed line) accretion. The "clumpy" gas accretion rate is shown by the black, dashed line. The present 
epoch, 2 = 0, occurs at 13.7 Gyr. 



is being delivered much closer to the galaxy disk (see also 
section [53111 ^^id hence appears to occupy a smaller spa- 
tial scale than the shocked gas. The gas that will later 
be accreted and shocked, on the other hand, has a much 
more uniform or spherical distribution. The shocked gas 
develops into a hot halo, much further from the disk. 

4. THE GROWTH OF THE STELLAR DISK 

Now that we have identified how gas gets to the galaxy, 
and whether it shocks or remains unshocked as it ap- 
proaches the disk, we wish to take it one step further 
and discover what role the various gas accretion mecha- 
nisms play in building the galaxy stellar disk. 

To identify disk particles, the galaxies are first aligned 
so that the disk angular momentum vector lies along the 
2-axis. ^z/ ^circ is calculated for each star particle in the 
galaxy, where Jz is the angular momentum in the x-y 
plane, and icirc is the momentum that a particle would 
hav e in a circular orbi t with the same orbital energy (see 
also lAbadi et al.ll2003[ ). Disk stars are identified as those 
having near circular orbits, so that ^z/icirc > 0.8. For a 
full discussion on our disk identification, see Brook et al. 
(2008, in prep). 

Fig. [S] shows the fraction of stellar disk mass at red- 
shift zero that has formed from gas of various sources, 
for each of our four main galaxies. Again, the total halo 
mass of the four galaxies is listed below their respec- 
tive bar, for reference. The total of the blue and red 
bars shows the fraction of stellar mass in the disk that 
formed from smoothly accreted gas, with red indicating 
the fraction that shocked as it entered the galaxy, and 
blue indicating the fraction that remained unshocked. 



The green bar shows the fraction of the stellar disk mass 
that formed from "clumpy" accretion, i.e., gas that at 
some point belonged to another galaxy halo. Finally, the 
black bar shows the fraction of stellar disk mass that was 
accreted directly as stars from other galaxies, while the 
white shows the fraction of stellar disk mass that formed 
from gas that was already in the galaxy at z = 6. 

Fig. [1] showed that smoothly accreted gas dominates 
over clumpy accretion for all four of our galaxies. Fig. [5] 
shows that this smooth gas goes on to form stars, so that 
the majority of the disk is composed of stars formed from 
smoothly accreted gas rather than from gas accreted in 
mergers or stars accreted in mergers. Table [3] is similar 
to Table[2]but for the stellar disk at redshift zero. Unlike 
Table [21 in which the fraction of smooth and clumpy gas 
totaled to f.O, in Table [3] the fractions of stars formed 
from smooth, clumpy, and "early" gas now sum to 1.0, 
where "early" gas is that which was already in the galaxy 
at z = 6, i.e., as far back as we trace the history of our 
gas particles. However, the fractions of stars that formed 
from shocked or unshocked gas (the last two columns in 
Fig. [5]) again sum to 1.0 (representing the total of stars 
formed from smoothly accreted gas). 

It can now be seen by comparing Tables [3] and [2] that 
the amount of shocked gas that makes it to the disk to 
form stars is lower than the total amount of shocked gas 
that is accreted to the galaxy. In fact. Table [3] shows that 
for all but our most massive galaxy, the stellar disk mass 
is dominated by stars formed from unshocked, smoothly 
accreted gas. This domination of unshocked gas in build- 
ing the stellar disk is even true for MWl, despite the fact 
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Fig. 4. — This figure shows histograms of the maximum temperature reached for all of the smoothly accreted gas for each of our four 
simulations (solid line). The dotted lines show the component of smoothly accreted gas that remains unshocked after being accreted to the 
galaxy, while the dashed line shows the gas that is shocked after entering the galaxy. 



TABLE 3 
Disk Star Forming Properties 




Fig. 5. Images from our MWl simulation at rodshifts 4 and 
3, identifying gas that will be smoothly accreted to the galaxy as 
either unshocked or shocked gas. Left column: The distribution of 
particles that will be (or have been) accreted but remain unshocked 
are marked in green. Right column: Particles that will be (or 
have been) shocked as they are accreted are shown in green. The 
underlying colors represent a gas density map, for reference, with 
black being least dense and white being the densest structures. 
The frames are centered on the main MWl progenitor at each of 
the redshifts shown. Each frame is ~1 comoving Mpc on a side. 
Faint, green boxes indicate Kvir at each time. 



that the majority of smoothly accreted gas to MWl gets 
shocked. 



galaxy 


^smooth 


^clumpy 


fshock 


funsliock 


H579 


0.96 


0.03 


0.00 


1.00 


DWFl 


0.75 


0.14 


0.06 


0.94 


MWl 


0.64 


0.30 


0.24 


0.76 


GALl 


0.63 


0.29 


0.69 


0.31 


H277 


0.67 


0.23 


0.32 


0.68 


MWl.lo 


0.70 


0.28 


0.28 


0.72 



Note. — The first two columns, famooth and 
clumpy 1 are the fractional quantities with re- 



spect to the total mass in stars in the galaxy 
disk at 2 = 0. Here, fdumpy includes only stars 
that formed from gas accreted from other ha- 
los (and not stars accreted directly as stars from 
satellites). By definition, famooth i fdumpy i and 
fcarly (not listed here) sum to 1.0. The last two 
columns, fghock and funahocki are the fractional 
quantities of stars that formed from shocked gas 
or unshocked gas with respect to the total stellar 
disk mass that formed from smoothly accreted 
gas that is in the galaxy at z = (and sum to 
1.0). 

4.1. The Early Growth of the Stellar Disk 

The importance of unshocked gas in the growth of the 
stellar disk is shown more explicitly in Fig. [7l which 
shows the star formation histories of those stars that are 
in the galaxy disk at redshift zero, for each of the four 
main galaxies (note, not the SFH of the entire galaxy). 
Here we can see that the disk SFR is, in fact, dominated 
at basically all times by stars forming from smoothly 
accreted, unshocked gas, for all but the most massive 
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Fig. 6. — Cumulative fraction of 2 = stellar disk mass by source, for each of our four galaxies, red: fraction of stellar disk mass 
that formed from shocked gas; blue: fraction of stellar disk mass that formed from unshocked, smoothly accreted gas particles that never 
belonged to a satellite galaxy before joining the main galaxy; green: fraction of the stellar disk mass that formed from "clumpy" accretion, 
i.e., gas that at some point was part of another galaxy halo; black: fraction of stellar disk mass that was accreted directly as stars from 
other halos; white: the fraction of stellar disk mass that formed from gas already in the main galaxy at 2 = 6. The total halo mass (in 
M0) of each galaxy is listed below its respective bar. 

the galaxy grows in time and reaches even higher virial 
temperatures, the radius at which shocking occurs in- 
creases (discussed further in Section [5331 . This increase 
in shock radius will potentially make later accretion to 
the disk even more difhcult, due to even longer cooling 
times. 

We have quantified the different delay times (the time 
between accretion to the virial radius and star formation 
in the disk) of the two types of smooth gas accretion. We 
find those stars that form from shocked and unshocked 
gas, and measure the average time between a gas parti- 
cle's accretion to the virial radius and the time when it 
begins to form stars. This delay is shorter for the un- 
shocked gas. The typical delay time for unshocked gas 
in our MWl and GALl galaxies is 1-2 Gyr. For shocked 
gas, the delay times in MWl are typically 3-4 Gyr, while 
for the hotter GALl the times are longer, 4-6 Gyr. 

Note that these delays are NOT intended to be repre- 
sentative of the cooling times of the particles. The delay 
between accretion and star formation includes both the 
cooling time of the gas particle onto the disk, the natural 
time the gas might spend in the disk before producing 
stars, and subsequent delays in SF due to feedback pro- 
cesses in the disk, which vary as a function of galaxy mass 
(Brooks e t al. 2007). In the two lowest mass galaxies, the 
delay time for unshocked gas is typically longer than for 
our two highest mass galaxies (3-4 Gyr rather than 1-2 
Gyr) due to the increasing importance of feedback effects 
at lower masses. 



galaxy, GALl. Even in GALl, unshocked gas dominated 
the growth of the stellar disk at an early period (to z ~1); 
only later does shocked gas begin to dominate. 

It is interesting to compare Fig. [7] to Fig. [3l which 
shows the gas accretion rate onto the virial radius with 
time. Note that in no case is the gas accretion rate ex- 
actly balanced by the star formation rate, though again 
we emphasize that we are only including disk SF here 
(and excluding, for example, the SFR of the bulge). Our 
two most massive galaxies are particularly interesting. 
Fig. [3] shows that their gas accretion histories are gener- 
ally dominated by smoothly accreted, shocked gas. How- 
ever, the star formation history of the MWl disk is in- 
stead dominated by smoothly accreted, unshocked gas. 
Star formation from gas that has been shocked initially 
increases steadily with time, until about 5 Gyr, after 
which the disk SFR due to this shocked component still 
remains relatively low, and eventually declines. GALl is 
a dramatic example of this. Shocked gas dominates the 
gas accretion history of the galaxy at essentially all times 
in the traced history. Yet again the initial star formation 
of the disk is instead dominated by unshocked gas accre- 
tion, while the star formation from shocked gas slowly 
increases with time, before peaking at about 9 Gyr and 
again declining. 

This delay in star formation for the shocked gas is not 
surprising. As this gas reaches higher temperatures (and 
entropies) , it has much longer cooling times, slowing the 
rate at which it can cool to the disk and form stars. As 
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Fig. 7. — The star formation history for those stars that are in the disk of each galaxy at z = 0. The star formation rates are broken 
down according to the accretion mode of the gas progenitor particles. Stars that formed from "clumpy" accreted gas are represented by 
the black, dashed line. The star formation rate of stars formed from smoothly accreted, unshocked gas is shown by the blue, dot-dashed 
line. Finally, the star formation rate for stars formed from smoothly accreted, shocked gas is shown in red (solid line). 2 = occurs at an 
age of 13.7 Gyr. 

The main conclusion, however, is that we expect longer 
cooling times for shocked gas compared to clumpy or 
unshocked gas. Thus, clumpy and unshocked gas will be 
able to make it to the disk to form stars at a faster rate 
than the shocked gas. In the next two sections we discuss 
the impact of cold flow (unshocked) gas accretion on disk 
formation. 



4.2. Observable Properties 

Fig. [5] shows the rest frame B band surface bright- 
ness for GALl aX z = 1. This image was created using 
Sunrise, a Monte Carlo radiative transfer code ( JonssonI 
l2006f ). The growth of the disk in GALl is dominated by 
unshocked gas accretion prior to z = 1, which clearly re- 
sults in an observable disk by this redshift. In a future 
paper we will investigate the evolution of disk sizes since 
z = 1. 

Fig. [5] shows Sunrise images of the / band surface 
brightness of stars at z = that form from the vari- 
ous accretion modes in our MWl simulation. The left 
panel shows stars formed from unshocked, smoothly ac- 
creted gas; the middle panel shows stars that formed 
from shocked, smoothly accreted gas; the right panel 
shows stars that formed in the galaxy from gas that had 
been accreted from other halos (i.e., clumpy gas accre- 
tion) . These panels show that the stars formed from the 
clumpy and unshocked, cold gas are found in both the 
stellar disk and bulge component, while the stars that 
are able to form from shocked gas are instead predomi- 
nantly in a disk. This disk formed from shocked gas is 
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Fig. 8. — The rest frame B band surface brightness images of 
GALl at 2 = 1, from Sunrise (in the case with dust scattering 
turned off). Note the obvious disk, even at this redshift. 



also thinner than the disk created by the unshocked or 
clumpy gas. The scale height of the shocked disk is ~25% 
smaller than the scale height of the unshocked disk. 

The morphology indicated by Fig. [5] can be explained 
in terms of the accretion times of the various gas com- 
ponents, and the subsequent cooling rate of that gas. 
The earliest disk stars are predominantly formed by un- 
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Fig. 9. — / band surface brightness images from Sunrise, for our MWl galaxy at z = (again with dust scattering turned off). Left: 
stars formed from unshocked, smoothly accreted gas. Middle: stars formed from shocked, smoothly accreted gas. Right: stars formed from 
clumpy accretion gas. 

formed from shocked gas have a fonger disk scale length 
(4.1 kpc) than those formed from unshocked (3.5 kpc) or 
clumpy gas (3.3 kpc). Although the picture is confused 
slightly due to the longer cooling times of the shocked 
gas, this scale length difference may not be surprising 
given the differing accretion times for these two compo- 
nents. Although there is still a low level of unshocked 
gas accretion at the virial radius all the way to redshift 
zero, the majority of the unshocked gas is accreted prior 
to z = 2 when the mean angular momentum is expected 
to be lower, and when chaotic merging will heat the disk. 
After this time, shocked gas dominates the gas accretion 
to the halo, when the mean angular momentum is ex- 
pected to be higher, and will therefore form a disk with 
a longer scale length. The angular momentum content of 
this gas will be investigated in detail in a future paper. 

4.3. Comparison with SAMs 

An important point to draw from the previous two 
sections is that, particularly as seen in our two most 
massive galaxies, a significant amount of star formation 
is possible in disks prior to z = 1 (when the universe 
was 6 Gyr old). As we have shown above, this early star 
formation is predominantly due to the contribution from 
unshocked gas accretion. Of course, the disks continue 
to grow since z = 1, but an existing disk is in place 
already by that time. This is in accord with observations 
that find that la rge galaxy disks i nust be assembled 
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Fig. 10. — Average surface brightness versus radius for the Sun- 
rise images in Fig. |9] Blue diamonds show the contribution from 
smoothly accreted, unshocked gas; red X's from smoothly accreted, 
shocked gas; green +'s from clumpy accretion gas; black triangles 
show the total. It can be seen that unshocked gas dominates the 
disk magnitudes at all radii, but that the shocked gas has a longer 
scale length (listed as h in the upper right corner). See the text for 
further discussion. 

shocked and "clumpy" gas, which dominate the accre- 
tion rate to the galaxy at the earliest times. This gas 
can rapidly cool and form stars, but at these early times 
the galaxy is still undergoing a series of mergers (major 
and minor) that will both heat the forming disk (lead- 
ing to a larger scale height for the stars created from 
unshocked and clumpy gas) and lead to bulge formation. 
The shocked gas, on the other hand, is both the dominant 
mode of accretion at later times and takes substantially 
longer to cool to the disk after accretion. The shocked gas 
will therefore settle into the disk after the chaotic merger 
period. This late settling allows the shocked gas to cool 
to a th in disk before forming stars (see also Bro ok et al.l 
l2004bf l. Hence, the scale height of the stellar compo- 
nent formed from this shocked gas is smaller than the 
unshocked case. 

Fig. [10] dissects these Sunrise images further, show- 
ing the radially averaged surface brightness profiles for 
each of these three components of gas accretion. As ex- 
pected from Fig. [71 the stars that form from unshocked 
gas dominate the surface brightness (shown as blue di- 
amonds; clumpy as green -|-'s; shocked as red X's; total 
as black triangles) . These profiles are fit well by a Sersic 
bulge component plus exponential disk. 

The scale lengths, h, for each type of accretion are 
listed in the upper right corner of Fig. 1101 The stars 
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Theoretical models of disk formation suggest instead 
that large disks must have formed since z = 1 and hence 
undergone a la r ge change in scale length since then, 
(|Mo et all [19981 : iMao et al.l[l99l Ivan den BoschI [19981. 
Even the most up-to-date model s still suggest a d iscrep- 
ancy with observational results (jSomerville et al.|[2008l ). 
with observed disks being larger than predicted at z = 
1. We will address disk sizes in a future paper, as it also 
requires a detailed examination of the angular momen- 
tum content of the accreted gas. We focus here on the 
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age of the stellar disk. 

It is particularly interesting to examine the stellar disk 
growth for our two most massive galaxies, which have 
been dominated by shocked gas accretion over most of 
the history we trace here. The unshocked gas accreted 
to these galaxies never reaches a temperature above 3/8 
Tyir- SAMs adopt a model in which this gas would first 
have been heated to the virial temperature of the halo 
as it entered the virial radius, thus delaying its ability 
to cool quickly to the disk. As demonstrated in Fig. 01 
the temperature difference between the shocked and un- 
shocked gas can be a factor of 10. This temperature 
discrepancy may lead to significant differences between 
the SFRs found here and in SAMs. 

In order to quantify the difference in SFRs between 
SAMs and our simulations, we estimate the effect on the 
disk SFH if the cold, unshocked gas had instead first been 
heated to T^r . The resulting change in the disk star for- 
mation rate is not straightforward to assess. Even our 
cold flow gas does not immediately form stars when it 
reaches the disk. Rather, there is a delay in star for- 
m ation due to feedba ck effects in the disk. As found 
bv lBrooks et al.l (|2007D . this delay is a function of galaxy 
mass, being longer in low mass galaxies, so that low mass 
galaxies form stars inefficiently compared to higher mass 
galaxies. 

As discussed in Section 14.11 the time between accre- 
tion to the virial radius and formation of disk stars is 
much longer for shocked gas than for unshocked gas. We 
find the shocked and unshocked gas particles accreted 
between each output time step interval, and the mean 
difference in formation time for stars that subsequently 
form from that gas. To quantify the affect of heating our 
unshocked gas to high temperatures, we adopt this av- 
erage difference in formation time between shocked and 
unshocked star particles, and add this delay to the for- 
mation times of the star particles born from unshocked 
gas. Although not precise, this method captures the rel- 
evant physics, including cooling times and suppression of 
star formation by feedback, and will provide a meaning- 
ful estimate of the effects of assuming that all disk gas 
has been shock heated. 

The result of adding this delay is shown in Fig. [Til 
for MWl and GALl. The top panels show the actual 
disk SFRs from our simulations as a solid line, exactly 
as in Fig. [3 but now showing the total disk SFR from 
unshocked, shocked, and clumpy accreted gas together. 
The dashed line shows what the resulting SFR would 
be if the unshocked gas had been heated to the same 
temperatures as the shocked gas. The middle shows the 
resulting cumulative mass of the stellar disk as a function 
of time, again as actually seen in the simulations (solid 
line) and with the false delay added to the unshocked gas 
(dashed line). The bottom panels show the difference 
between the two disk masses in the middle panel (the 
delayed disk mass as a fraction of the actual disk mass 
from the simulation). 

Fig. [Tl] demonstrates that if the cold gas had reached 
virial temperatures that the SFR in the disk would have 
been significantly delayed. Rather than being able to 
build the stellar disk within the first few Gyr, the disk 
does not begin to grow substantially until z ~1. The 
bottom panel shows that the disk is 50% less massive 
until z = 1 in the case of strong heating than it actually 



is in the simulation, where the unshocked gas does not 
approach the virial temperature. 

5. NUMERICAL ISSUES 

We discuss in this section the impact of key numerical 
effects on our results. 

5.1. Time Resolution 

The effect likely to have the largest impact on our re- 
sults is the interval between our outputs. Given our finite 
step interval of ^--^320 Myr, it is possible that a particle 
with a short cooling time may radiate away its energy 
quickly enough that we miss an increase in temperature 
and entropy associated with a shock event (making our 
shocked fractions a lower limit). This omission is more 
likely to occur for those particles shocking near the disk 
(at ^ 30 kpc) where the densities are much higher (and 
hence cooling times shorter) than at the virial radius. 

To investigate the effect of our step size on our results, 
we have run a new simulation with output steps at 1/12 
of the time resolution presented above. This new simu- 
lation, II277, is a Milky- Way like galaxy with a rather 
quiet merger history and a total halo mass of 7.1x10^^ 
Mq (see Table [ij at z = 0. For the same step inter- 
val as was used above, the gas accretion results for this 
galaxy agree well with the trends shown for the previous 
four galaxies. From Table [2] it can be seen that the frac- 
tion of shock accreted gas lies nicely between DWFl and 
GALl, as might be expected for this mass of galaxy. 

A change in time resolution has an effect at two differ- 
ent stages of our analysis. First, as just mentioned, we 
may miss particles with short cooling times that should 
have been counted by our "shocked" definition. Second, 
"clumpy" gas is identified as belonging at any output 
time step to a halo other than the main galaxy halo. We 
are likely to miss gas particles that belong to another 
galaxy for less than '^320 Myr (the time between output 
steps), e.g., because they are subsequently lost from the 
galaxy due to SNe. The result of missing this clumpy 
accretion will be to increase the number of smoothly ac- 
creted particles instead. 

These two stages of our analysis are not indepen- 
dent of one another, since the sample of shocked gas is 
drawn only from those gas particles that are identified as 
smooth accretion. We wish to quantify them separately. 
We first test how the number of particles identified as 
"clumpy" changes as a function of time resolution. We 
investigate the change at 12x the original step interval 
used above (i.e., fine sampHng corresponds to ~27 Myr 
resolution per output). 

When using our standard step interval, the smooth 
gas makes up 73% of all the gas accreted by galaxy II277 
over the history of the trace. By increasing the output 
interval resolution by a factor of 12, the gas identified 
as smoothly accreted declines, so that it composes 66% 
of all the gas accretion over the galaxy's history, for a 
total decline of 7%. Thus, by definition, the fraction of 
gas accreted as "clumpy" over the galaxy's history must 
increase by 7%. 

This result means that our absolute numbers stated 
above for the fraction of gas that is smoothly accreted 
to each galaxy is an upper limit. An examination of 
the accretion history of this galaxy shows smooth accre- 
tion with characteristic spikes indicative of mergers. The 
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Fig. 11. — The top panels show the disk SFR in our simulations, as in Fig. [7] but now summed together as the solid line. The dashed 
line shows the effect on the disk SFR if unshocked gas took as long to reach the disk and form stars as the shocked gas. MWl is shown 
on the left, and GALl on the right. The middle panels show the cumulative stellar disk mass in each case. The bottom panels show the 
difference of the two disk masses shown in the middle panel, i.e., the delayed disk mass as a fraction of the disk mass from the simulation. 



amount of gas in these spikes identified as "smooth" de- 
creases in our fine time step limit, but the spikes still re- 
main, suggesting we are still missing some gas in mergers. 
By generously estimating the mass of gas contributing to 
these accretion spikes, we find that we are missing only 
an additional 1% of clumpy gas accretion Thus , we agree 
with previous studies fKOS. lMurah et al.ll2002l ) that find 
that the majority of gas mass is accreted smoothly rather 
than in mergers. 

Next, we use the "smooth" particles identified in the 
standard output resolution (320 Myr outputs, where 
smooth gas accounts for 73% of all the gas accreted over 
the history of the galaxy) to isolate the effect of output 
interval on the fraction of particles identified as shocked. 
At finer step intervals, the fraction of smooth gas that is 
shocked increases by 2.3%. This increase is quite small, 
and occurs primarily at 2 > 2, likely to due to the fact 
that we do not adequately sample cooling times at high 
redshift. To verify this, we calculated the instantaneous 
cooling rates for particles in an annulus between 30-40 
kpc (as close to the disk as we trace them) and inside 10 
kpc of Kvir at each of our "standard" analysis intervals. 
For H277, the virial radius does not reach 30 kpc until 
z = 2.65, and the cooling times of the particles within 
Kyir are lower than our standard step resolution of 320 
Myr until this redshift. After this redshift, however, the 
galaxy grows quickly and the cooling times are longer 
than our output intervals, so that we are capable of eas- 
ily identifying shocked gas. Thus, it is not surprising 
that the addition of shocked gas when increasing time 
resolution comes at high redshifts. 



We have done a similar cooling time analysis for our 
four main galaxies. In the case of DWFl, MWl, and 
GALl, we do not fully resolve cooling times within the 
halos until redshifts 3, 3, and 4.5, respectively. However, 
because DWFl is not massive enough to support a stable 
shock until z ^0.6, our time step undersampling at high 
z has no impact on our resulting fraction of shocked gas. 
For MWl and GALl, we can expect about a 5% change, 
as we found in h277, and perhaps less since we resolve 
their cooling times to higher redshift than in the H277 
case. 

As an interesting note, our lowest mass galaxy, H579, 
always has cooling times longer than our standard out- 
put step interval (though it is never massive enough to 
shock, so this is not important for our results). It may 
seem odd at first that this galaxy has such long cooling 
times, but it is a result of the lower densities in this halo. 
Our UV background is capable of ionizing much of the 
hydrogen gas in the halo of this galaxy, causing it to be 
inefficient at cooling. We verified that turning off our UV 
background lowers the instantaneous cooling times sub- 
stant i ally, by almo st an order of magnitude (|Efstathioul 
fl99l iQuinn et al.l fT9961 . 

5.2. Entropy Generation 

By definition, a shock is a phenomenon that leads to 
a sharp discontinuity in the physical properties of a gas. 
Due to the nature of smoothed particle hydrodynamic 
codes, these properties for any given gas particle will be a 
smoothed estimate based on the properties of its nearest 
neighbors. Thus, there is some concern that SPH will do 
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a poor job of resolving shocks. 

Gasoline ado pts the standard SPH artificial viscos- 
ity formulation (jMonaghanl Il992f ). with a correction 
to det ect and reduc e viscous transport iri shearing re- 
gions (Balsara"l995D. IWadslev et all (|2004f ) investigated 
Gas oline ' s perf ormance in this area with both the stan- 
dard ISodI (|1978D sho ck tube test, and a spherical adi- 
abatic collapse test ()Evrar dl fl98l . It was found that 
Gasoline does an excellent job of matching the ana- 
lytic solutions to within a few particle separations, and 
matches the post-shock entropy values very well. 

It can be possible for low numerical resolution to arti- 
ficially broaden shocks (due to entropy generated by the 
artificial viscosity used in SPH), allowing gas to radiate 
energy while passing through a shock rat her than heating 
and undergoing an entropy jump ( Hutchings fc Thomad 
|200C1() . This broadening is unlikely to be a n issue at the 
high resolutions used here (|SDringe]|l2005D . but we ver- 
ified numerical convergence by using a lower resolution 
version of our MWl run, MWl.lo, and performing the 
same analysis. This run has 1/8 of the mass resolution 
and 1 /2 the force resolutio n of MWl (and is the same 
galaxy discussed in detail in lGovernato et ani2007t l. Ta- 
ble [2] demonstrates that the shocked fraction in this run 
is in excellent agreement with MWl. Since our shock 
definition includes a temperature criterion, this implies 
that we are not missing high temperature gas, and thus 
that we have a sufficient numerical resolution to identify 
entropy jumps and are not missing heating in shocks. 
Hence, our cold flow gas accretion is not an artifact of 
numerical resolution. 

5.3. An Alternate Shock Definition 

Previous groups who have attempted to investigate the 
role of cold flows in galaxy gas accretion have used a dif- 
ferent criterion to identify shocked and unshocked gas, 
adopting a strict temperature cutoff. We examine in this 
section how our results differ if we use the same temper- 
ature definition as other authors. 

K05 and OPT08 have noted that the temperature dis- 
tribution of accreted particles can be bimodal, with a 
minimum at 2.5x10^ K. They use this minimum to sep- 
arate gas into cold accretion and hot accretion modes. 
We also see this minimum in the temperature distribu- 
tion, as shown in Fig. [12] for our MWl galaxy. This 
figure shows the maximum temperature of gas particles 
that have been accreted to the galaxy between successive 
redshift intervals. Up until z = 0.5, there is a strong 
bimodality. The dashed vertical line shows the sepa- 
ration at 2.5 xlO^ K used by K05, K08, and OPT08. 
Our accreted gas particles clearly demonstrate this same 
minimum. These two groups have previously noted that 
this minimum is a function of the cooling curve, and has 
therefore has little direct relation to the actual modes of 
hot or cold accretion. However, it is a convenient point 
at which to separate the gas into hot or cold accretion 
modes. 

To assess how our results are affected by these defini- 
tions, we repeat the above analysis to identify shocked 
gas in our smoothly accreted gas component, using the 
temperature criterion of K05 and OPT08. The same 
set of smoothly accreted particles is used, so that the 
only change between the two definitions will be in the 
resulting shocked or unshocked fraction (or, more appro- 



priately, hot and cold fraction). The difference resulting 
from these two separate definitions is remarkably small. 
The temperature criterion leads to slightly more hot gas 
being accreted to the galaxy than shocked gas, but the 
largest difference is only 7%, in the case of GALl. The 
remaining galaxies have differences of 3% or less. H579 
never reaches the temperature cutoff identified by K05, 
and thus both studies agree than it has never accreted 
shocked/hot gas. 

These tests lead us to conclude that the two differ- 
ent definitions yield very similar results. This agreement 
may be partly due to the fact that we also impose some 
sort of temperature criterion (requiring that our gas not 
only undergo an entropy increase, but reach 3/8 T^^^ as 
well). However, this result was not guaranteed. We em- 
phasize that the temperature criterion alone is based on 
the physics of the cooling curve, and does not reflect the 
physics of the gas accretion shocks. Hence, we choose to 
add an entropy criterion to our deflnition, as it is unique 
to each galaxy halo and the physics of the various accre- 
tion modes. 

5.4. The Role of the UV Background 

iGovernato et al.l (|2007t l found that the mass function 
of our galaxies is resolved down to a minimum of 64 DM 
particles. For the runs used here, that corresponds to to- 
tal halo masses of 8.1x10^ Mq, 6.1x10^ Mq, 4.9x10"^ 
Mq, and 1.9x10*^ M© for H579, DWFl, MWl, and 
GALl, respectively. These masses are well below the 
limit where the UV fiel d wil l suppress gala xy formation 
(IGovernato et al.l l2007t jEfeta thiou 199^ lOuinn et al " 
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20061 ). and thus we avoid issues due to sub-resolution 
merging. That is, there should be no unresolved halos 
that arc being counted incorrectly as "smooth" accre- 
tion. This assertion was verified by running the same 
analysis on a lower resolution version of our MWl galaxy, 
MWl.lo, with a mass resolution of 1.9x10® Mq at 64 
DM particles. This mass is just below the limit that the 
UV field should suppress galaxy formation. Indeed, the 
smooth gas fraction agrees to within a few percent with 
MWl (see Tabled]). MWl.lo does have a five percent 
higher fraction of smooth accretion, but given fluctua- 
tions in the simulations, this is excellent agreement. 

The dominant effect of the UV background in our sim- 
ulations is to suppress galaxy formation at low masses, 
and hence lower the amount of clumpy accretion. Gas 
which would have otherwise been in small halos is now 
free to be smoothly accreted to our galaxies. This was 
confirmed by running MWl.lo with the UV background 
off (and cooling and SF turned off as well). The results 
for this galaxy, MWl. ad, are shown in Table|2] It is clear 
from the first two columns that the amount of "clumpy" 
accretion material increases from about 20% to about 
50%. This confirms that our UV background is partly 
responsible for reducing the amount of gas accreted from 
other halos, and hence increasing the amount of gas avail- 
able for smooth accretion. 

The more interesting effect of the UV background is on 
the amount of shocked gas accretion to the galaxy. Our 
UV background (in combination with the mean baryonic 
density of the universe as a function of redshift) sets a 
temperature fioor in the IGM, and hence an entropy fioor 
for IGM gas. For a galaxy to be capable of shocking ac- 
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Fig. 12. — This figure shows histograms of all of the smoothly accreted gas between successive redshift intervals for our MWl galaxy. 
The dashed vertical line shows the temperature divide used to separate hot and cold accretion modes used by K05 and OPT08. 



creted gas, the virial entropy of the galaxy halo must 
be above this background entropy by a certain factor. 
Hence, based on our definitions, the UV background sets 
a minimum halo mass that will be able to shock, as func- 
tion of time. 

5.5. Cooling 
5.5.1. The Adiabatic Case 

iBirnboim fc Dekell ()2003f ) and DB06 showed that the 
shock radius in the case of adiabatic accretion (i.e., no 
cooling) is well matched to the virial radius in a spheri- 
cally symmetric case. They also demonstrated that when 
cooling is turned on, the same galaxy is not initially mas- 
sive enough to support a stable shock, so that accreted 
gas shocks at the disk of the assembling galaxy. Only 
after a critical mass is reached can the galaxy support a 
shock. 

We have again used MWl. ad with the UV background, 
cooling, and SF turned off (i.e., the adiabatic case) to test 
these results. As noted above, the amount of smooth 
accretion drops for MWl. ad, as much of the gas is able 
to remain bound to its dark matter halo. However, from 
Table [2] (last two columns) it is clear that the amount of 
smooth gas accretion that is shocked also increases in this 
case, increasing to 84%. This behavior is expected, and 
is similar to the standard spherical model in which all gas 
gets heated to the virial temperature of the galaxy as it 
is accreted. The fact that there is unshocked gas at all 
is due to the fact that this simulation is not spherically 
symmetric, and filaments of cold gas still exist that are 
capable of penetrating well inside of the halo without 
shocking. 




Fig. 13. — These contours show the radii as which particles 
are identified as shocked, for our MWl.lo and adiabatic, MWl. ad, 
runs. Black represents the lowest concentration of particles at a 
given time, and white peaks are the highest. The solid lines show 
the virial radius as a function of time. 

Fig. [T3] shows the radii at which particles are shocked 
inside of Kvir, as a function of time, for the MWl.lo and 
adiabatic, MWl. ad, runs. Contours are used for clarity, 
with black representing the least dense concentration of 
particles, and white peaks the most dense. The solid 
lines show the virial radii of the galaxies as a function 
of time. The main difference between the two runs is at 
high redshift. In the adiabatic case, particles are able 
to shock inside the virial radius even at high z, while in 
the MWl.lo run there is no shocking at the highest red- 
shifts. Shocks begin to occur near the disk at redshift 4, 
and shocked particles fill the space between the disk and 
the virial radius at lower rcdshifts. The radius with the 
highest concentration of shocked particles moves outward 
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with time in both cases. 

5.5.2. Metal Line Cooling 

Metal line cooling has recently been added to Gaso- 
line (Shen et al., in prep.). MWl.lo was rerun with 
metal cooling on, to find if the new cooling has any ef- 
fect on the results here. However, the smoothly accreted 
gas (be it shocked or unshocked) has nearly primordial 
metallicities, and the resulting fractions of shocked and 
unshocked accreted gas are identical to the results for 
MWl.lo with only H and He cooling. Likewise, the frac- 
tion of disk stars formed from shocked and unshocked 
gas is also identical to MWl.lo. 

6. DISCUSSION & CONCLUSIONS 

We have tracked the history of gas accreted to four high 
resolution disk galaxies that span two orders of magni- 
tude in galaxy mass. For all of these galaxies, early disk 
growth can be predominantly attributed to cold gas, and 
up to Milky Way masses the SF of the stellar disk is 
dominated by unshocked, cold flow gas accretion at es- 
sentially all times. This dominance of cold gas in the 
formation of disk stars can exist even when the the ma- 
jority of gas accreting to the virial radius is being shocked 
and reaching high temperatures. At galaxy masses above 
L*, an early phase of disk growth prior to z = 1 can be 
attributed to the presence of cold gas accretion. 

We have examined gas temperature histories in rela- 
tion to the virial temperature of the halo that they are 
accreted to. While the shocked component is capable of 
reaching temperatures near the virial temperature of the 
halo, the unshocked component of gas accretion remains 
below 3/8 of the virial temperature at all times. The 
important role of unshocked, cold flow gas in the growth 
of galaxies is consistent with previous studies that adopt 
an explicit temperature cutoff, determined by the cool- 
ing curve, to identify cold flow gas. We confirm that 
the ability of smoothly accreted gas to shock is a strong 
function of mass. However, our focus here has been to 
study the subsequent role of this cold gas accretion in the 
growth of stellar disks. Because the unshocked gas does 
not reach the same high temperatures as shocked gas, it 
is capable of cooling to the disk faster and forming stars. 

It is important to note that this early cooling of un- 
shocked gas to the disk is not a manifestation of the 
historic "overcooling" problem. This can be seen by 
examining the results from our stellar mass - metal- 
licity relationship, which matches the observed rela- 
tion at both low and high redshifts (iBrooks et al.ll2007t 
iMaiohno e t al.''2008l . As discussed bv both lBrooks et al.l 
(|2007^ and Maiolino et al. (2008), this match is achieved 
by overcoming the overcooling problem with both high 
resolution and a physically motivated SN feedback 
model. Previous works consistently produced metallic- 
ities too high when compared to observations at high z, 
due to overcooling and a rapid consumption of gas and 
subsequent metal enrichment. Instead, our SN feedback 
scheme regulates the SFRs of our galaxies, preventing 
them from too quickly consuming their gas. The early 
SF seen in our more massive galaxies is due instead to 
the accretion of cold gas that can cool rapidly to the disk 
to form stars. 

Despite running our simulations in a fully cosmologi- 
cal context that includes galaxy mergers, the dominant 



gas supply is from smoothly accreted gas that has never 
belonged to another galaxy halo. For a Milky Way mass 
galaxy with a quiescent merging history like the galaxies 
shown here, only ~25% of the accreted gas is acquired 
from other galaxy halos, resulting in a similar fraction of 
the disk stars forming from this gas. 

It is our goal to assess the impact that cold flow 
gas accretion may have in altering current models of 
disk galaxy formation. A direct comparison with semi- 
analytic models is difficult to do and beyond the scope of 
this paper. Current SAMs include the possibility of rapid 
gas cooling for low mass galaxies. It has been suggested 
that this regime is s imilar to the cold fl o w gas accretion 
seen in simulations (|Croton et al.ll2006l ). ICattaneo et al.l 
(2007) per formed a detailed comparison of the GallCS 
SAM (Ha tton et al.t 2003) to the SPH resuhs of K05, and 
found that the cold gas accretion rates in the two schemes 
were in agreement. This implies that the amount of cold 
gas currently available in SAMs may agree well with sim- 
ulations, though this needs to be confirmed with further 
comparisons. 

However, the results in this paper may emphasize a 
separate problem within the cold flow regime. Even in 
the case of rapid cooling in SAMs, the temperature of the 
gas is initially set to the virial temperature of the galaxy 
halo, and cooling times are based on this assumption. 
Prior to the development of a hot halo, the free fall time 
determines when the gas reaches the disk in either the 
SAMs or the simulations. However, SAMs do not yet 
model the impact of fllaments on galaxy growth. Our 
two most massive galaxies, which are ab le to develop a 
hot h alo at high redshifts (4 to 6; see also lDe Rossi et all 
l2008l ). are still capable of accreting cold gas that never 
approaches the virial temperature of the halo in fllaments 
at high z. These fllaments allow a signiflcant amount of 
cold gas mass to accumulate inside of the galaxy, even 
for massive galaxies capable of supporting a shock at 
early times. This cold gas cools quickly to the disk to 
form stars. This rapid settling to the disk allows for the 
building of stellar disks at higher redshifts than predicted 
if all gas is heated to the virial temperature, as in the 
standard model. 

We have shown that the discrepancy in temperatures 
between the shocked and unshocked gas components 
leads to signiflcantly longer cooling times for the shocked 
gas (note that their densities will be discrepant as well, 
since most of the unshocked component are found in 
dense filaments, exacerbating the difference in cooling 
times) . If the unshocked gas had reached the virial tem- 
perature, the star formation rates in the disk would have 
been delayed (in the cases shown here, until after z = 
1). A study of the angular momentum content of the 
accreted gas will be examined in future work, but the 
presence of cold flow gas accretion, particularly in flla- 
ments in massive halos, allows for the building of large, 
massive disks at high z, in agreement with observations. 
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